Critical theory of the topological quantum phase transition in the AKLT— SZH chain 
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We systematically study the phase diagram of S = 2 spin chain by means of density-matrix 
renormalization group and exact diagonalization. We confirm the presence of a dimer phase in 
the AKLT-SZH model and find that the whole phase boundary between dimer and SZH phases, 
including the multicritical point, is a critical line with central charge c = 5/2. Finally, we propose 
and confirm that this line corresponds to SO(5)i Wess-Zumino-Witten conformal field theory. 
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In recent years, investigations of topological phases 
and phase transitions has attracted great attention in 
condensed matter physics* and quantum information 
theory^. Topological phases are characterized by a bulk 
gap separating excitations from the ground state and 
by the presence of gapless edge modes. A topological 
phase cannot be deformed continuously into a conven- 
tional, topologically trivial phase without going through 
a phase transition, in which the gap closes and the edge 
modes merge with the bulk. The quantum Hall state^ 
and the recently discovered topological insulators^bre ex- 
amples of topological states of quantum matter. 

Topological phases even appear in one dimensional 
(ID) systems. The Haldane phase^ of integer quantum 
spin chains is an example of a symmetry protected topolog- 
ical phas^®. In 1987, Affleck, Kennedy, Lieb, and Tasaki 
(AKLT)P introduced a family of exactly solvable integer 
spin models with valence bond solid (VBS) ground state 
and proved that VBS states share the key features of 
Haldane gap liquids for integer spin chains. The ground 
state of AKLT model can be exactly formulated in terms 
of Schwinger bosonsP, |1< AKLT ) = f[ {ij} (a\b] - b\a^) s \0), 
where 5 is the local spin. The VBS states are unique 
ground states of the AKLT Hamiltonian for periodic 
boundary condition (PBC); for open boundary condi- 
tions (OBC), however, the ground state is (5 + l) 2 -fold 
degenerate with edge spins 5/2. The parent Hamilto- 
nian for the VBS states |\f rAKLT ), here for 5 = 2 is 
conveniently defined in terms of projection operators^, 

^AKLT = J2K 3 P 3 (i,j) + K A P A (i,j), (1) 

(ij) 

where K 3 , K4 > and P3 and P4 project onto the spin 3 
and spin 4 subspaces at the bond respectively. The 

edge spin is boson-like 5 = 1 and the square of the time 
reversal operator T satisfies T 2 = 1. 

In 1998, Scalapino, Zhang, and Hanke (SZH) intro- 
duced an SO(5) symmetric superspin modePwith an ex- 
act VBS ground state. The local spin in the SZH model 
transform under the five-dimensional vector representa- 
tion of SO (5). SZH presented an exact ground state wave 
function expressed as a matrix product state of Dirac T 



matrices, and showed that the ground state is 16-fold de- 
generate for OBC. This model can be naturally mapped 
onto a spin 2 chain, with fermion-like 5 = 3/2 edge spinal 
and T 2 = — 1. Following SZH's work, more VBS states 
with higher symmetry groups have been constructe d 11 ! 12 !. 

The ground state of SZH model is given by |<3/ SZH ) = 
X; {mi} Tr(r mi r m2 ...r mN )|mi---m N ) , where the T m 

fulfill r a r b = 2S ab + 2iT ab and m, is a vector label of 
the SO(5) group, which can also be interpreted as the 
rrii = —2,-1,0,1,2 quantum numbers of spin 2. The 
corresponding Hamiltonian is given by 

tfSZH = ^ j 2 p 2 ( i;i ) + J A P 4 (i,j), (2) 

(ij) 

where J 2 , J\ > and P 2 and Pa again are the SU(2) bond 
projection operators onto spin 2 and 4, respectively. The 
ground state is unchanged up to an SO (5) rotation. 

Like for topological insulators, the bulk topology is 
related to the edge states of an open chain. Conse- 
quently, AKLT and SZH models describe different topo- 
logical phases. As recently pointed outPI, odd integer 
spin AKLT models are protected by a couple of symme- 
tries, e.g. time reversal, while even integer spin AKLT 
models are solely protected by global SU(2) symmetry. 
Since the first case is characterized by half-integer edge 
spin, the SZH model is similar to odd spin AKLT models 
(notice that the spin 3 AKLT model exhibits an edge spin 
3/2 like SZH). Both, the AKLT and SZH phases are thus 
protected by symmetries, but the SZH phase seems to be 
much more robust. Given the topological distinction of 
the two ground states, we construct a model Hamiltonian 
interpolating between the AKLT and SZH models: 

H(a) = (1 - a)P AKLT + aP SZH (3) 
= ^[aP 2 (i,i) + (l-a)P 3 (i,j)+^4(i,i) • 

(ij) 

Here, we set J 2 = K 3 = 1, J 4 = K 4 = (3. As the edge 
state is robust unless the gap closes, there must exist one 
or several topological quantum phase transitions (TQPT) 
when varying a from to 1. 
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FIG. 1: (color online) Ground state phase diagram of the 
Hamiltonian (3} obtained by DMRG with N = 600 sites. The 
location of second order phase transition is indicated by the 
red line. The region above the red line belongs to SZH phase, 
while the region below the dashed blue line belongs to AKLT 
phase. The intermediate regime between the topologically 
different SZH and AKLT phases are dimer phases. 



In Ref.[TU the ground state phase diagram of the 
AKLT-SZH model (^J has been studied by means of 
density-matrix renormalization group (DMRGJPl Be- 
tween the AKLT and SZH phases, a possible dimer phase 
has been found. The most interesting feature is the pres- 
ence of a multicritical (MC) point, at which a direct 
TQPT occurs between the two distinct phases. How- 
ever, there are still some important questions. (1) Does 
the dimer phase indeed exist in the thermodynamic limit? 
(2) What is the phase diagram for f3 < 0.27? (3) What 
is the central charge and low-energy effective theory at 
the MC point? (4) Does the effective theory spread over 
the whole phase boundary? 

Motivated by these questions, in this paper we will 
revisit the model ^ and try to answer all open questions. 
For the present study, we keep up to m = 3400 DMRG 
states with more than 16 sweeps to get converged results, 
and the truncation error is less than 10 -9 . We use both 
OBC and PBC with system sizes up to N = 600 sites. 

Dimer order parameter. We confirm the presence of 
the dimer phase for /? > 0.27, see the phase diagram 
Fig.[T] In addition, we find a much larger dimer phase 
for (3 < 0.27. To clarify the presence of the dimer phase, 
we have calculated the dimer order parameter (DOP), 
V = |(SjSj+j.) — (Sj+iSi+a)! i as a function of a and /3, 
where i = N/2 and we imposed OBC. In Fig.ETa), ex- 
emplarily the DOP as a function of a at (3 = 0.2 for 
different system sizes is shown. In the inset, the finite- 
size scaling is also performed to obtain the DOP in the 
thermodynamic limit. As shown in Figs.[T] and [2] when 
j3 decreases from /? = 1, also the maximum amplitude 
of DOP decreases, and the dimer phase shrinks with a 
monotonously. Finally, at the MC point f3 = 0.27, the 
DOP vanishes for any a. Surprisingly, for smaller val- 
ues of j3 there is a revival of the dimer phase, which be- 




FIG. 2: (color online) Dimer order parameters (DOP) of the 
Hamiltonian (J3J with different of a and /3. (a) DOP for /3 = 
0.2 at system size N = 40, 60 and 100, as well as in the 
thermodynamic limit N = oo. The finite-size scaling is given 
in the inset, (b) DOP for different f3, which disappears at the 
multi-critical point ft — 0.27. 

comes even larger and the magnitude of DOP increases, 
indicating a very robust dimer phase at small j3 re- 
gion. In the limit f3 — 0, the P4 term vanishes and 
we expect a huge ground state degeneracy. Therefore, 
it will be extremely difficult to get converged DMRG 
results and we only consider /3 > 0.1 in the phase di- 
agram Fig.[T] Like for the general spin 1 chain, we ex- 
pect a generic dimer phase containing the Hamiltonian 
.ffdimer = J2i ~Po{h * + 1) with an enhanced SU(5) sym- 
metry. We performed DMRG calculations to verify this 
guess. We also studied the two Hamiltonians interpo- 
lating between -ffdimcr and the dimer phase for j3 > 0.27 
and for j3 < 0.27, respectively. In both cases, we found in 
the second derivative of the ground state energy with re- 
spect to the interpolation parameter a signal for a second 
order phase transition. So it might be that the general 
spin 2 chain contains three different dimer phases. The 
nature of the phase transition between dimer and AKLT 
phases seems to be a higher order (i.e., third or more) 
phase transition. Therefore, it will be extremely difficult 
to locate the exact phase boundary numerically. 

String order parameter. To distinguish different 
topological phases, one way is to use the hidden 
string order parameter (SOP) defined as Q(6) = 
(Aiexp(i6J2 j kZ!iAk)Aj). For spin-S 1 AKLT model, the 
peak of SOP is at C = n/S for A\ = Sf. But this 
SOP cannot be used to distinguish the AKLT and SZH 
phase since it is finite in both phases. Fortunately, Tu 
et o/P proposed a SOP with A, = i(5f )[(5f) 2 - 1] and 
9 C = n which knows about SO(5) symmetry of the SZH 
model and which is finite only in the SZH phase. There- 
fore, we can use the SOP to characterize the TQPT be- 
tween SZH and other phases (in addition to d 2 E/da 2 ). 
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FIG. 4: The 100 lowest levels of the entanglement spectrum 
— log Xi for system size N = 100 at /3 — 0.2. Xi are the 
eigenvalues of the reduced density matrix. 



FIG. 3: (color online) String order parameters of H(a) for (a) 
/3 = 0.2 and (b) /3 = 1.0, for system size N = 40, 60 and 100, 
and the extrapolated value for N = oo. 



In Fig.|3j the SOP for /3 = 0.2 and /3 = 1.0 is shown for 
several system sizes and for N — > oo. In both the AKLT 
and dimer phase, G(tt) — > 0. In the SZH phase, however, 
G(~k) remains finite, and approaches its maximum value^ 
4/5 2 = 0.16 in the SZH limit a = 1. 

Entanglement spectrum. To identify different phases, 
we can also use the entanglement spectrum (ES), which 
was originally introduced by Haldane^ in the context of 
FQHE and later applied to spin chains in a momentum 
basis^. The ES is defined as the set of eigenvalues of the 
reduced density matrix pa = TrB|^')(\E'|, with A being a 
subsystem and B the remainder of the system. Here we 
use the real space ES with PBC to characterize the dif- 
ferent phases. As shown in Fig.[4j the lowest level of ES 
in AKLT phase is 9-fold degenerate, this degeneracy is 
associated with the edge spin 1. By increasing a, a finite 
size splitting of the 9 ground states occurs (chain length 
N = 100). After entering the dimer phase, the lowest 
value of the ES becomes non-degenerate. By further in- 
creasing a, the system enters the SZH phase, the lowest 
level of the ES becomes degenerate again, but now the 
degeneracy is 16-fold in agreement with edge spin 3/2. In 
all the three phases, the lowest level of ES is separated 
by a large gap from the other levels. 

Central charge. The most interesting feature is the 
presence of the MC point, as shown in Fig.[T] To char- 
acterize the critical theory we should determine cen- 
tral charge and scaling dimensions. The central charge 
can be obtained by calculating the von Neumann en- 
tropy of a subsystem A with length x, defined as 
5 vN = — Tr(pAhi/9A)- For critical systems, it has been 
established^ that 5 vN = (c/3) ln(x') + c x for PBC, and 
S vn = (c/6)ln(2a;')4-ln(g)+ci/2for OBC, where c is the 
central charge of the conformal field theory (CFT), c\ is a 
model dependent constant, and g is Affleck and Ludwig's 
universal boundary ternn^. For finite chains, we can use 
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FIG. 5: (color online) (a) Finite-size scaling of the peak posi- 
tion a c (N) of d 2 E/da 2 as a function of a at /3 = 0.27. (b) Von 
Neumann entropy S vN at the multi-critical point for different 
chain length with the fitted central charge c. Inset: S vN for 
different parameter points is shown for comparison. 



the conformal mapping x — > x' = (N/tt) sm(irx / N) . Us- 
ing this formula the central charge can be extrac ted in 
excellent agreement with the CFT predicitorPSHH. 

By performing finite-size scaling of the peak position 
of d 2 E/da 2 as a function of (3 and a, we find the exact 
position of the MC point at {(3 c ,a c ) = (0.27,0.8518). In 
Fig.[5ja}, the finite-size scaling is shown for j3 = 0.27. 
In Fig.p3[b), we show S vN at the MC point for different 
system sizes. The regression fit for 5' vN shows very good 
convergence with system size, indicating a central charge 
c = 5/2. For comparison, in the inset of Fig.[5]jb), we 
show the results for some points away from the MC point 
where S' vN starts to saturate indicating the opening of a 
gap. In the same way, we find c = 5/2 for the whole 
critical line between SZH and dimer phases. 

Effective field theory. To find the corresponding CFT, 
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TABLE I: List of CFT with central charge c = 5/2. S2/S1 is 
the ratio of the second and first primary scaling dimensions. 





SO(5)i 


SU(2)w 


SU(2) 2 x U(l) 


SU(2) 4 x Ising 


MC 


s 2 

Si 


1.6 


2.67 


2.67 


2.67 


1.56 



we have listed all the simple corresponding candidates 
with central charge c = 5/2 in Table]!] S2 / Si is the ratio 
between the scaling dimensions of the second and first 
non-trivial primary fields. To characterize the CFT the- 
ory at the critical line, we rescale and match the lowest 
three finite-size energy levels obtained numerically by ex- 
act diagonalization (ED) for systems with up to L = 12 
sites to the form of the spectrum of a CFT^U, 

E(L)=E 1 L+^(~ + h + h) . (4) 

Here the velocity v is an overall non-universal scale fac- 
tor and the scaling dimensions h + h take the form 
h = h + n,h = h + n, with n and n non-negative 
integers. h° and h° are the holomorphic and antiholo- 
morphic conformal weights of primary fields in the CFT. 
The momenta (in units 2tt/L) are such that k — h — h 
or k = h — h + L/2. We consider the energy spectrum at 
k = with scaling dimension S — 2h (i.e., h = h). We 
find that S2/S1 = 1.56 at the critical line. Comparing 
this value with the values in Table[IJ we find that only 
SO(5)i is compatible with our data, while all the other 
candidates can be explicitly ruled out. To further con- 
firm this, we also calculate the ratio (1 + S2) / (1 + Si) at 



k = 2it/L. We find (1 + S 2 )/(l + Si) = 1.27, which is 
also consistent with the value (1 + 1)/(1 + 5/8) = 1.23 
(here n = 1) of SO(5)i CFT, by considering the finite- 
size effect. We conclude that the critical line is described 
by SO (5)i Wess-Zumino-Witten (WZW) model. 

This result is particularly interesting since we expect 
three integrable models in the phase diagram of the 
general spin 2 chairP^ the Takhtajan-Babujian chain 
(SU(2) 4 WZW with c = 2), the permutation operator 
(SU(5)i WZW with c = 4), and the Reshetikhin model. 
From the latter model it is believed that it could be de- 
scribed by SO(5)i WZW with c = 5/;Pl It turns out 
that this model is located near our critical line in the 
phase diagram of the general spin 2 chain. Indeed we 
find c = 5/2 for the Reshetikhin model, but we leave its 
connection to our critical line as an open question^. 

In conclusion, we have confirmed the phase diagram 
of the AKLT-SZH chain; remarkably, also for (3 < 0.27 
a dimer phase is present. The critical line separating 
the SZH and dimer phases has the same central charge 
c = 5/2 and can be described by SO(5)i WZW theory. 
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